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Abstract 

We discuss theoretical aspects of parton distribution functions for very high energy 
scattering in relation with upcoming measurements in DIS and hadron-hadron colli- 
sions. 

1 Introduction 

The Large Hadron Collider will operate with very high parton luminosities. As these 
rise steeply for decreasing momentum fraction x, a large number of events sample 
gluon and sea-quark distributions at x <C 1. Understanding the theoretical accuracy in 
the determination of these distributions is relevant both to obtain reliable predictions 
for cross sections of hard processes and to investigate new aspects of QCD physics at 
very small x. 

Current determinations of parton distribution functions for x < largely de- 
pend on deep inelastic scattering data. We start the discussion in Sec. 2 by comment- 
ing on such determinations, with a view in particular to the forthcoming longitudinal 
structure function measurements. 

Most of the phenomenological pdf analyses place cuts on the low-Q 2 region, in 
order to ensure the applicability of perturbation theory. However, a lot of the exper- 
imental information onx C 1 physics lies at present with data that involve scales of 
few GeV. In Sec. 3 we turn to discuss how s-channel methods, designed to describe 
high-energy scattering down to lower and lower Q 2 , might be used in the context of 
the parton picture. Measurements of jet production in the DIS current fragmentation 
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and central region could be exploited in this framework to extract information on the 
nonperturbative parameters of multiple parton scattering. 

We illustrate in Sec. 4 the use of the above methods for estimating power cor- 
rections to structure function's evolution from multiple scattering. We emphasize 
the distinctive behavior of the results in the region of intermediate scales just above 
1 GeV. Through momentum sum rules and evolution, the low-x and low-Q 2 region 
may affect predictions for processes at much higher x and higher mass scales. Some 
additional comments on related issues are given in Sec. 5. We conclude in Sec. 6. 



2 Flavor- singlet parton evolution 

Present determinations of parton distribution functions are based on fits to available 
collider data [H El El H] , using renormalization-group evolution equations to connect 
measurements at different mass scales \i. The region x 1 is dominated by flavor- 
singlet evolution, 

d ( fs \ _ ( Pqq Pqg\ ^ ( fs 



dln^ 2 V fg J V P m Pgg ) ® \ fg ) ' ^ 

where f g and f$ are the gluon and sea-quark distributions, and is the pertur- 
batively calculable evolution kernel. The f g distribution in the range x < 10 -2 , in 
particular, is relevant for measurements of production processes dominated by gluon 
fusion at the LHC, including heavy flavor production El [7] and Higgs boson pro- 
duction 0[8]. 

At present the extraction of / 9 at x < 1CT 2 depends on DIS data for the Q 2 
derivative of the structure function F 2 (F 2 = dF 2 /dha.Q 2 ), 

F 2 ~f s , F 2 ~ P qg ® f g [1 + 0(A 2 /Q 2 )] + quark term . (2) 

The kernel used for standard pdf determinations is evaluated in fixed order of 
perturbation theory (NLO and, in some of the analyses now becoming available [HE], 
NNLO). As indicated by Eq. (j2J), however, the extraction of f g is especially sensitive 
to higher- loop corrections to the gluon — ► quark splitting kernel P qg . In particular, for 
x < 1 logarithmically enhanced contributions a k s +l x~ - l ln^ 1 x to P qg are present 
for any k > 1, and need to be resummed. The numerical impact of terms of this type 
on the global fits is examined e.g. in [I]. 

While the summation of In x terms in the pure-gluon sector has been studied ex- 
tensively for quite some time (see reviews in 0, [10], and references therein), only 
recently have the first analyses appeared [HI [TJl H3] that implement both gluon [H] 
and quark [T5] next-to-leading lnx corrections to the matrix kernel Py. This is the- 
oretically appealing, as quark corrections are required in order to merge consistently 
the x — > expansion with the short-distance behavior and the renormalization group, 
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and opens the possibility for resummed analyses to become directly relevant to phe- 
nomenology 

The production of 6-quarks at the LHC will be sensitive to f g at x ^ 1CT 2 . 
Recall that the theoretical uncertainties on NLO predictions for the b cross sections, 
while on the order of a few ten percent at the Tevatron [16] . increase to well over 
a factor of 2 at LHC energies 13 [16]. Improving upon present predictions will 
likely involve a variety of physical effects, from higher-loop corrections beyond NLO 
to nonperturbative processes both in the initial and the final states. An improved 
treatment of x — > contributions for both the pdfs and the short-distance cross 
sections should help understand the sources of the large uncertainties at the LHC, 
and possibly reduce them. 

The production of Higgs bosons may also receive sizeable contributions from low-x 
gluon events at the LHC, depending on the Higgs mass range and the kinematical 
region in the Higgs rapidity and px- See e.g. [T7] for studies of high-energy effects on 
the accuracy of Higgs boson cross sections. 

The xCl region will be probed experimentally by the forthcoming measurements 
of the longitudinal structure function Fl [18J. A recent phenomenological study of 
Fl may be found in Ref. [12]. This presents a thorough comparison of fixed-order 
predictions through NNLO [19] and predictions including the lnx resummation for 
the Fl coefficient functions [15]. The analysis p2] shows that i) the NLO and NNLO 
corrections to Fl are large, as observed in [19], and lead to strong instabilities causing 
perturbation theory to be badly behaved at low x; ii) resummed results improve this 
behavior, and allow one to obtain more reliable theoretical predictions throughout a 
wider range in the kinematical variables. This is potentially significant, since the Fl 
measurements will provide an independent observable, to be combined with F 2 , to 
perform a complete flavor-singlet analysis and probe the accuracy of the theory at 
x < 10 -2 more stringently than ever so far. 

A different, but possibly related issue concerns the initial conditions at low mass 
scales for pdf evolution. Working to NNLO, Ref. [20] observes that rather different 
features than in [T2"j \W\ |2"T] are obtained in the fixed-order predictions for Fl if 
different assumptions are made [201 E2] for the initial conditions. We note that this 
marked dependence on the initial condition may be related to using a fixed-order 
truncation of perturbation theory that is not well-behaved at low x, and could likely 
be reduced in the improved theory including resummation. 

On the other hand, while resummed evolution schemes provide better theoretical 
control on the x <C 1 region [TTl [T21 [TB"] . note that the coefficients themselves of 
resummed perturbation theory signal the onset of dynamics beyond the a s expansion. 
One way to see this is to examine the singularity structure of resummed coefficients in 
the plane of the effective x — ► anomalous dimension 7^(0^) [23]. The strong branch- 
point singularities in gluon-channel coefficients [23] i? u (a s ) ~ l/y/l/2 — 7 a ,(a s ), with 
7^ — > 1/2 for x — > 0, imply that effects beyond the perturbation expansion are to be 
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included in order to fully describe the high-energy limit (see also discussion in [24]). 

In Eq. d2J) we have indicated explicitly contributions to F 2 suppressed by powers 
of l/Q 2 - Because present DIS data that are relevant to extract the gluon density 
for x < 1CT 2 do not have very high Q 2 , these power corrections may be expected 
to have non-negligible effects on the estimate of the theoretical accuracy on pdf's 
for the LHC In the next sections we turn to methods for the subset of power-like 
contributions that comes from graphs with multiple gluon scatterings. 



3 Relating parton approach and s-channel approach 

Methods to take account of multiple scatterings are based on the s-channel picture 
of deep inelastic collisions. See [25] for an introduction to the physical picture, 
and [261 E3 EHl EH] for motivation from high-density QCD. The main advantage 
of this approach is that it gives a formulation that can be used down to small Q 2 , 
incorporating nonperturbative physics at very high energies in Wilson-line operator 
matrix elements. 

The s-channel approach leads to a different picture of the hard collision than the 
parton approach, as it works in a different reference frame, and uses different degrees 
of freedom. Nevertheless, the two pictures are not incompatible. In the region where 
their domains of validity overlap, they must describe the same physics. Ref. [30] 
presents a framework to connect the two pictures. 

The method [30j is based on constructing explicitly an s-channel representation for 
the renormalized parton distribution function in terms of Wilson-line matrix elements. 
In this representation the quark distribution f q is given by the coordinate-space con- 
volution (Fig. [[)) 



xf q (x,fj,) — J dz J db u(fi, z) E(z,b) — UV , (3) 
where H is the hadronic matrix element of eikonal-line operators, 
E(z,b) = J[dP'] (P'\± Tr{l-V\b + z/2)V(b- z/2)}\P) , 

{r+oo -\ (4) 

-ig J dz-A+{0,z-,z)t a [ , 

z is the transverse separation between the eikonal lines, b is the impact parameter, 
and the function «(//, z) is evaluated explicitly in [30] at one loop using the MS scheme 
for the renormalization of the ultraviolet divergences z —>■ 0. 

The convolution ([3]) expresses the fact that for x 1 the operator defining the 
quark distribution function creates the partonic system, made of a color-triplet eikonal 
line and an antiquark, at large longitudinal distances y~ ~ 1/ (xP + ) far outside the 
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Figure 1: Quark distribution function in the s-channel picture. 



target [25]. The MS result [30J can be recast in a physically more transparent form 
in terms of a cut-off on the z integration region (see also [31j ) , as long as the scale \i 
is sufficiently large compared to the inverse hadron radius: 

xf q / P (x, fj) = ^ \ dh % 9 ( z2 ^ > fl2 ) b ) ' ( 5 ) 

where a is a renormalization scheme dependent coefficient [30] . 

a = 2e 1/6 " 7 « 1.32657 , (6) 

with 7 the Euler constant. The Wilson-line matrix element b) receives contribu- 
tion from both long distances and short distances. At small z it can be treated by 
a short distance expansion. By using renormalization-group evolution equations, the 
leading term in the expansion can be related to a well-defined integral of the gluon dis- 
tribution function [30]. At large distances z, S has to be fit to data, or parameterized 
by models consistently with bounds from unitarity and saturation J29J [36] . 

Besides parton distribution functions, the s-channel representation discussed above 
is potentially interesting for jet production. To this end, we note that Eqs. 
imply introducing a cut-off in rapidity [SU] in order to separate slow and fast par- 
tons in the picture of Fig. [TJ and to factor out the Wilson-line matrix element. The 
dependence on the method of carrying out this separation enters in higher orders 
for the quark distribution, and in leading order for the gluon. Both the quark and 
the gluon channel contribute to the jet structure of the final states. Nevertheless, if 
one considers jet cross sections by fixing e.g. the total minus momentum of the dijet 
system, the sensitivity to the rapidity cut-off is suppressed, as long as the jets are 
produced sufficiently far from the fragmentation region of the target. Then the ap- 
plication of the s-channel method to measurements of jet production will give rather 
direct information on the contribution of multiple parton scatterings. 



Note that current shower Monte-Carlo generators include model parameterizations 
of multiple interactions in order to produce realistic event simulations, see e.g. [32| [33]. 
In this respect, jet rapidity distributions in the DIS current fragmentation and central 
regions, once analyzed using the s-channel picture, are potentially useful in order to 
constrain the nonperturbative parameters of such models. 

The physical picture described in this section relies on the large separation in 
lightcone coordinate between creation of the partonic system and interactions with 
the target. This separation is of order a hundred fermi in the proton rest frame 
for x < 10~ 2 (see remark after Eq. (j4j)). The condition for the applicability of the 
approach is satisfied in the case of high-energy collisions of large nuclear targets as 
well. The approach should be relevant for the physics of high-p t probes in heavy ion 
collisions at the LHC and of nuclear parton distributions |34j . 



4 Power corrections from the s-channel 

The representation ([3]), evaluated in a well-prescribed renormalization scheme, is the 
key ingredient that allows one to relate [301 ES] results of s-channel calculations for 
structure functions [26j [29], [36] to the OPE factorization, 

F 2 = C®f + ^ 2 C^®f^ + ... , (7) 

and, in particular, to identify power-suppressed contributions to the Q 2 evolution of 
F 2 of the form 

1 A 2 

F 2 ~ P gg ®f g [1 + 5} + quark term , 5 ~ ^ a k {a s — —) k . (8) 

k>i X ^ 

The enhanced x — > behavior in the power correction 8 in Eq. ([8]) is produced 
from graphs with multiple gluon scatterings, and is consistent with observations of 
approximate geometric scaling in low-x data [37J. 

Refs. [351 [38] introduce moments A 2 of the matrix element H in Eq. ([4]), 

^Hj^/^^r 1 JdbE(z,b) , (9) 

and express the l/(Q 2 ) h correction in terms of A 2 (fc) times coefficients computable as 
functions of a s , x and \nQ 2 , schematically in the form 

m _( dF 2 \ ~ A 2 (A;) 



dlnQ 2 \d\nQ 2 J LP ^ (Q 2 ) 
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The first term in the right hand side is the leading-power parton result, and the 
moments A in the subleading terms are dimensionful nonperturbative parameters, to 
be determined from comparison with experimental data. Observe that any number 
of rescatterings contribute to the moment A 2 (fc) through the eikonal operators (jlj). 
Thus the correction of order l/{Q 2 ) k receives contribution from the exchange of ar- 
bitrarily many gluons. This is characteristic of the small-x power correction and can 
be contrasted with counting rules arguments, valid at large x, that link the power in 
Q 2 with the number of parton lines exchanged in the t-channel. 

The result of determining the nonperturbative A parameters from F2 data [39] 
at both low and high Q 2 is shown in Fig. [2] in the left hand side plot [35] . The 
corresponding power correction is plotted on the right hand side of Fig. [2j Here the 
correction is normalized to the full answer and multiplied by —1. 
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Figure 2: (left) The result of fitting the A 2 parameters to the data[39]; (right) power 
corrections to dF 2 /d In Q 2 versus Q 2 at different values of x|35j. 

The analysis [33] indicates that with physically natural choices of the parameters 
in the nonperturbative matrix elements in Eqs. (jlJ),(jSD one can achieve a sensible 
description of data for x < 1CT 2 in a wide range of Q 2 and still have moderate power 
corrections to dFi/d In Q 2 . Corrections turn out to be negative and below 20 % for 
x > 10~ 4 and Q 2 ^ 1 GeV 2 . This observation suggests that the power expansion 
is not breaking down, and should still work at least to the values of x considered in 

Fig. m 

However, we also see from Fig. [2] that for small x the corrections fall off slowly in 
the region of medium Q 2 , Q 2 ~ 1 — 10 GeV 2 , behaving effectively like 1/Q U with v 
close to 1 [35]. For instance, one has v ~ 1.2 for the curve x = 10~ 3 in the right-hand 
side plot of Fig. [2j As a consequence of the slowly decreasing behavior, the power 
corrections stay on the order of 10% up to Q 2 of a few GeV 2 for x < 10 -3 . This 
slow fall-off differs from parameterizations of higher twist commonly used in global 
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analyses (see e.g. Ref. and may be relevant for phenomenology as it affects the 
medium Q 2 region of the data that are useful to extract f g at low x. 

5 Further comments 

The region of low x and low-to-medium Q 2 influences predictions at much higher x 
and mass scales through momentum sum rules and evolution. We give below some 
further comments. 

It is worth noting that the biggest contribution to the power correction to F2 
in [35] comes from the longitudinal component. The derivative dF^jd In Q 2 has 
smaller power corrections than dF 2 /d In Q 2 . This provides additional motivation for 
the separate measurement of F^ |18j . Similar conclusions are reached from a different 
perspective based on the fits [H [21], investigating the effect of power-like terms on 
global analyses at both high x and low x. (These fits also include parameterizations 
of power effects from self-energy graph insertions. We recall here that only relatively 
few results are known as yet on such effects in flavor-singlet observables. See [40J for 
a recent study, and early discussions in [IT].) 

The curves in Fig. [2] are obtained using NLO parton distributions. It is natural 
to expect a change in the power correction when going from NLO to NNLO. We note 
that the decrease in the low-x gluon at NNLO [3] is consistent with the possibility 
that NNLO partons give smaller power corrections. 

Observe that in the calculation described in Sec. 4 the slow fall-off with Q 2 results 
from summing the terms proportional to the moments A in Eq. ([9]). These in turn 
are obtained from expanding the s-channel answer in powers of 1/Q 2 , and enforcing 
consistency with the standard parton framework. It will be of interest to compare 
the Q 2 behavior found here with the behavior due to the anomalous dimensions that 
result from nonlinear evolution equations such as those, e.g., used in [261 123 [28j HQ]. 

It will also be of interest to investigate the relation of the result fl9l). f|T0l) for the 
power correction with the x-rescaling form proposed in [42] and applied to nuclear 
targets. This should likely involve the finite lightcone-time cut-off that enters in 
the high-energy eikonal approximation (see e.g. [43]). reflecting the fact that the 
eikonalized projectile-target interactions do not spread out to arbitrarily large times 
in the far past and the far future. 

The behavior observed in Fig. [2] suggests that the power corrections may be char- 
acterized by a nonperturbative scale substantially larger than Aqcd- It is possible 
that this can be related to the dynamical cut-off on large transverse distances z im- 
posed by unitarity requirements ("black disc" limit) on the correlator H in Eq. To 
pursue this, further studies of the moments are warranted. Also, the analogue of 
these moments for eikonal operators ([4]) in the adjoint representation will be relevant 
for processes coupled directly to the gluon distribution. An especially interesting case 
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is that of diffractive DIS events [44J , where cross sections depend quadratically on the 

5 correlator, and gluonic contributions typically dominate quark contributions by one 
order of magnitude. See comments in [38], and references therein, on the possible role 
of the "blackness" cut-off in hard- diffractive data. 

6 Conclusion 

The partonic structure of protons and nuclei will be probed for small values of x 
at the LHC. DIS measurements of longitudinal Fl, expected in the coming year, 
will provide very valuable, new experimental input. As discussed in Sec. 2, improved 
theoretical tools for the evolution of parton distributions are starting to become avail- 
able that include the results of next-to-leading In x resummation for both gluon and 
quark channels, as is required by consistency with the short-distance behavior and 
the renormalization group. 

Much of the present information on the region x < 10~ 2 comes from deep inelastic 
data involving scales of few GeV. Physics beyond the perturbation expansion is likely 
to play a role in this regime. The physical picture of hard collisions that is designed 
to go down to lower and lower Q 2 when x is small is the s-channel picture. A method 
to use results of s-channel calculations in the context of the parton framework is dis- 
cussed in Sec. 3, and applied in Sec. 4 to corrections to structure function's evolution 
suppressed by powers of l/Q 2 but enhanced as x — ► 0. 

The results presented build a physical picture of the sea-quark distribution for very 
high energies, that allows one to discuss quark saturation. Nonperturbative power- 
like effects are expressed in terms of moments of Wilson-line eikonal correlators. The 
picture suggests the possibility of using measurements of jet leptoproduction in the 
intermediate rapidity region to measure effects of parton rescatterings. 
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